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All nuclei with even numbers of protons and of neutrons have ground states with zero angular 
momentum. This is ascribed to the pairing force between nucleons, but simulations with random 
interactions suggest a much broader many-body phenomenon. In this Letter I project out random 
Hermitian matrices that have good quantum numbers and, computing the width of the Hamiltonian 
in subspaces, find ground states dominated by low quantum numbers, e.g. J = Q. Furthermore I 
find oAA-Z, odd-A'^ systems with isospin conservation have relatively fewer J = ground states. 

PACS numbers: 21.10.Hw,21.60.Cs,24.60.Lz 
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Most quantum many-body systems cannot be exactly 
solved, even numerically. Because generic many-body 
systems are complex and classically chaotic, one way 
to model the spectrum is through random matrices 
These matrices must be Hermitian, of course, but there 
are no other exact symmetries. 

Real many-body system often have non-trivial symme- 
tries such as rotational invariance and isospin invariance, 
which give rise to states with exact quantum numbers 
such as angular momentum J, J^, and isospin T, T^. 
But because the above random matrices do not have such 
symmetries, investigators tended to consider statistical 
properties of states with the same quantum numbers. 

In contrast, the ordering of different quantum numbers 
in spectra was associated with details of the interaction. 
For example, the fact that nuclei with an even number 
Z of protons and an even number N of neutrons always 
have ground states with angular momentum J = 0, while 
odd-Z, odd-iV nuclei frequently do not, was attributed to 
the pairing interaction [2, 0] • 

It was therefore a shock to discover that rotationally 
invariant but otherwise random two-body Hamiltonians 
tend to yield ground states with J = 0, just like 'realistic' 
interactions, even though such states are a small fraction 
of the total space j^-iS,]. This phenomenon is robust, in- 
senstive to details of the distribution of matrix elements 
0, occurs not only for fermions but also for bosons 0, 
and is relatively insensitive to the particle rank of the 
interaction Q. Over the past decade there have been 
many papers proposing explanations. As the distribution 
of many-body systems with two-body interactions tend 
to have a Gaussian distribution of stat es |lCll| . a number 
of authors have focused on widths Uj, , while others 
have statistically averaged in a single j-shell the coupling 
of multiple angular momenta [l3]. As a recent Letter 
stated, 'the simple question of symmetry and chaos asks 
for a simple answer which is still missing [9| . ' 

In this paper I return to random matrices and impose 
symmetries, first U(l) then SU(2). I show explicitly how 
combining 'internal' degrees of freedom with projection of 
good quantum numbers leads naturally to subspaces with 
small quantum numbers having the greatest widths, and 



thus dominating the ground state. This simple picture 
applies with equal ease to both fermionic and bosonic 
systems, and aside from subspace dimensions is indepen- 
dent of the detailed microphysics, helping to explain the 
robustness of the phenomenon. And considering two si- 
multaneous SU(2) symmetries, angular momentum and 
isospin, I find conservation of isospin in odd-Z, odd-A^ 
system changes dimensions of subspaces in such a way 
as to decrease the fraction of J = ground states, a 
prediction confirmed with detailed simulations. 

I start with U(l) symmetry and consider a wavefunc- 
tion which is periodic + 27r) = '0(</') which 
is an eigenstate of a general eigenvalue equation: 



(1) 



Without significant loss of generality I assume the wave- 
function -0 and the Hamiltonian 0') to be real. Her- 
miticity requires that H{(j),(f)') = H{ip',(j)), while U(l) 
invariance suggests that H can only depend on the dif- 
ference of angles: (/>') = F{(j) — (/)/). Combining Her- 
miticity with periodicity leads to F{x) — F{2ti — x). 

Inasmuch as is a periodic function, I make a Fourier 
decomposition, keeping in mind that F is real. 
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which can be inverted 



2 cos{mx)F{x)dx 
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SO that the integral is only over unique values of F{x). 

If F{x) is a randomly distributed variable about a; = 0, 
with a variance independent of x, then on average the 
value of hm is zero and the variance is easily computed: 



4 cos^(TOa;)F^ {x)dx 



4cos'^{mx)a^dx = 27r(l -\- Sm,o)'^^ 
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If the only degree of freedom is (j), then the hm are 
the eigenvalues of H , and m both labels the solutions 
and their symmetry. But now suppose there are (dis- 
crete) internal degrees of freedom, for instance if one has 
a many-body system that has an overall U(l) symmetry. 
In that case 'P{x) is a matrix-valued function of x, and 
h„i is a random symmetric matrix, with the dimensions 
of F counting internal degrees of freedom. 

If the matrix elements of F(a;) each have a variance of 
a independent of x, then the variance of the individual 
matrix element of are still given by ([4]). Hence the 
matrix for m = has twice the variance of m > and 
the ground state will be likely have m = 0. 

I can illuminate this by discretizing (j). Then the sym- 
metry forces a matrix of the form 



H 



( A B C D 
B A B C 
C B A B 



C B\ 
D C 
E D 



(5) 



\ B C D E ... B A J 



Adding 'internal' degrees of freedom means that 
A, B,C, . . . are now random symmetric matrices of the 
same dimension. Numerical calculations with ([5]) verify 
the accuracy of ([3]) and (j4]) and the dominance of m = 
quantum numbers for the ground state. 

Now consider SU(2), rotational invariance. Using the 
angles from spherical coordinates, the Hamiltonian is 
of the form H{9' (j)' ^Oif), but imposing rotational invari- 
ance means H can only depend on the angle u be- 
tween 6',(j)' and 0,0 as given by cos a; = cos f? cos 6*' -|- 
sin^sin^' cos((/) — 4>'). Then much like U(l) 



H{e'4>\ec^)=F{u:), 
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where F{uj) = F{uj + 2'k) is a periodic function and, using 
Hermiticity (and assuming H is real) F{uj) = F{2t: — uj) 
is symmetric with respect to w = tt. Expanding 



2J + 1 
47r 



Pj (cos w) 



(7) 



M=-l 



SO clearly the hj are again the eigenvalues, with Yjm as 
eigenfunctions and with the eigenvalues independent of 
iV, as one expects. 

Once more I assume F(w) to be a matrix-valued func- 
tion, and thus hj to be a symmetric matrix given by 

hj = 27r / Pj(cos ti;)F(w) sinwdw. (8) 

^0 

As before, let a be the variance of the matrix elements 
F independent of w. Then I estimate the variance of the 
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FIG. 1: (color online) Distribution of ground state quantum 
numbers for eight neutrons in the lp-0/ shell. Empty bars are 
the fraction of states in the model space with a given J, hor- 
izontally striped (red) bars are the fraction of ground states 
with a given J predicted by a random matrix (RM) model, 
and vertically striped (blue) bars are fraction of ground states 
with a given J from configuration-interaction (CI) diagonal- 
izations of an ensemble of random two-body Hamiltonian in 
a shell-model basis. 



matrix elements of hj 



a] = 47r2a-2 ( p] {cos uj) sin^ 
Jo 



(9) 



Eqn. ([9]) can be computed numerically, and leaving off 
the factor Air'^a'^, the values are 1.571, 0.393, 0.245, 0.178, 
0.139 for J = 0, 1, 2, 3, 4 respectively. 

This suggests that in many-body system, subspaces 
with low-valued quantum numbers will have larger 
widths. But in realistic, finite many-body calculations, 
subspaces with different Js have different dimensions. 
Furthermore, in the above argument each hj has in- 
dependent random matrix elements, which typically has 
a semi-circular density of states [ij, yet for many-body 
systems with only two-body interactions the density of 
states tends towards a Gaussian (lo|. 

I can approximately correct both deficiencies. First, 
following standard results on matrices with Gaussian- 
distributed matrix elements [1, I let 



^cff 
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be the width of the subspace of states with angular mo- 
mentum J. Then, for each J, I simply create Nj energies 
via a random Gaussian distribution of width tr}'^, and 
ultimately determine the fraction of ground states with 
angular momentum J. 

Finally, I compare against a variety of concrete sim- 
ulations via configuration-interaction calculations, that 
is, diagonalizing the Hamiltonian for fixed numbers of 
particles in finite single-particle spaces. Figure [1] shows 
the case of eight fermions (neutrons) in the lpi/2-lp3/2- 
O/5/2-O/7/2 shell- model space (which in nuclear physics 
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TABLE I: Ground state quantum numbers for single-species 
systems, comparing the percentage of ground states of a given 
J for configuration-interaction simulations (/ci) against the 
percentage computed using the simple random-matrix model 
(/rm) described in the text. As input, /apaco is the fraction of 
the states with a given J in the CI model space. The two cases 
are 8 fermions in a j = 21/2 shell, and the N = 7 interacting 
boson model (IBM). 
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corresponds to ^*Ca with an inert ^°Ca core), while Ta- 
ble U considers two more cases, eight identical fermions in 
a j = 21/2 shell, and six bosons in the interacting boson 
model (IBM) 14]. 

For each of these many-body systems I take an en- 
semble of rotationally invariant, two-body but otherwise 
random interactions (typically a few thousand cases), and 
tabulate the fraction fci of states that have a given an- 
gular momentum J. This should be compared with the 
native fraction of states with that J in each many-body 
space, /space = Nj/Ntot (^tot is the total dimension of 
the many-body space), and as noted originally the frac- 
tion with J = is dramatically enhanced. 

I also compare with the fraction of ground states with a 
given J predicted by my simple random matrix picture, 
/rm- The only input are the dimensions Nj and the 
universal variances computed in ^ and used in (jlOp . 

For such a simple picture, the random matrix model 
yields qualitatively excellent results, generally predicting 
the enhancement or suppression of different Js in the CI 
simulations relative to the native fractions /space- In par- 
ticular, not only does the RM model successfully predict 
an enormous enhancement of J = in the ground state, 
it predicts a mild enhancement of J = 2. 

This analysis suggests the predominance of angular- 
momentum zero ground states is primarily a function 
of the width of the angular-momentum-projected many- 
body Hamiltonian; furthermore, the width is largely de- 
coupled from the microphysics, instead depending only 
on the projection integrand (|9]) and on the dimensional- 
ity of subspaces with good quantum numbers. The sim- 
plicity and decoupling from the microphysics may be why 
the phenomenon is so robust and so universal. 

So far I have only considered angular momentum. Yet 
in nature, nuclei with even numbers of protons and even 
numbers of neutrons always have J — ground states 
while those with odd numbers of protons and odd num- 
bers of neutrons often do not, and if one runs ensembles 
of configuration-interaction simulations with two-body 
interactions that conserve angular momentum and and 



TABLE II: Distribution of ground state angular momentum 
J for systems with equal numbers of protons and neutrons, 
given as a percent for configuration-interaction simulations 
(/ci) against fraction in the simple random-matrix model 
(/rm). The columns marked 'conserved' means isospin conse- 
vation was enforced, while for those marked 'broken' isospin 
was maximally broken. The single-particle model spaces are 
Is-Od or (sd) and lp-0/ or (pf). 





conserved 


broken 


J 


J space 


J space 


/rm 


/ci 


/space 


/rm 


/ci 






{sd),Z -- 


= iV 


= 3 









0.8 


1.6 


42 


15 


3.6 


72 


32 


1 


2.5 


4.3 


32 


34 


10 


16 


31 


2 


3.5 


6.4 


17 


9 


15 


10 


14 


3 


4.2 


7.3 


7 


26 


17 


2 


15 


4 


4.1 


7.2 


0.5 


1.6 


16 


.1 


4 






(p.f),z-- 


= N 


= 3 









0.7 


1.2 


41 


11 


2.6 


66 


28 


1 


2.1 


3.4 


27 


36 


7.6 


16 


26 


2 


3.0 


5.2 


19 


7 


11 


13 


16 


3 


3.8 


6.2 


11 


23 


13 


3 


15 


4 


3.9 


6.6 


2 


2 


14 


0.5 


3 



isospin, this scenario is broadly reproduced: one gets a 
predominance of J = ground states for even-even cases 
but greatly reduced for odd-odd. 

I now generalize the above simple model, and consider 
widths that depend upon both total angular momentum 
J and total isospin T. If Njt is the dimension of a CI 
space with fixed J, T, then let the width be 

^fx = VnJ^^J^t (11) 

where ctt is also taken from 

Table |ll] shows the dimensions for several proton- 
neutron cases, both for fixed J, T and for fixed J alone, 
as well as the random matrix prediction for the fractions 
fjT and fj. For even-even cases (not shown to save 
space), there is little difference, but for odd-odd, break- 
ing isospin dramatically enhances the fraction of J = 
ground states. This can be considered a prediction of 
the simple random matrix model. The results can be 
traced directly to the subspace dimensions (here /space) - 
Specifically, as one goes from isospin breaking to isospin 
conserving, a smaller fraction of J = states go to T = 
than do J = 1, and the relative decrease of dimensional- 
ity of J = 0, T = makes the difference. 

I include the results of an ensemble of CI simulations, 
which verify the qualitative predictions of the random 
matrix model. While that quantitative agreement is sig- 
nificantly less than for the single-species case, the quali- 
tative trends do agree. 

In computing (jH]), I assumed each F(cj) to be uncor- 
related, so that the widths add incoherently. Suppose 
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FIG. 2; (color online) Fraction of ground states with angu- 
lar momentum J as a function of the correlation length (see 
EQ. I12p . for N — Z — 4 in the sd shell. Diamonds correspond 
to the TBRE result; only the J = 1 is shifted significantly, 
see dashed arrow. 



instead that that F(a;) and F{uj') are correlated, with a 
correlation length 7; then the width ([9]) becomes 

/>7r />7r 

47r^(T^ / dw I dw'Pj (cos w) sin ojPj (cos oj') sin w' 
Ja Jo 

xexp(-|tj-tj'|/7)-(12) 

In the limit 7 — )■ cx), which corresponds for a perfectly 
correlated, i.e. constant, F this vanishes for J > 0. 
Fig. ([2]) shows how the fraction of ground states with 
J varies with correlations length 7, for the case of ^''Mg. 

Even if F{uj) are not correlated, the hj still have corre- 
lations because / Pj(cosa;)Pj'(cosw) sin^ wdw ^ 0; this 
correlation was found via a much more complicated prior 
analysis [l^]- One immediate result follows from the in- 
tegral: while J = is correlated with all even J, it is 
not correlated with odd J. Further consequences will be 
pursued in future work. 

In summary, detailed quantum many-body simulations 
using ensembles of random interactions have shown sur- 
prising trends, most notably the predominance of angu- 
lar momentum J = in the ground state. To address 
this question, I have argued how one can project ran- 
dom matrices with good quantum numbers. Using only 
simple, universal integrals and the dimensionality of sub- 
spaces with good quantum numbers, I can qualitatively 



reproduce the features of ensembles of many-body sys- 
tems. In particular I find an dominance of J = for 
the ground state, although I reproduce other trends as 
well. When I further consider systems with two species 
(protons and neutrons) the random matrix model pre- 
dicts, and CI simulations confirm, that the subspace di- 
mensions of isospin-conserving systems suppresses J — 
ground states in odd-Z, odd-iV systems. 

The results of the model are qualitative, not quantita- 
tive. On the other hand, the effectiveness of the quali- 
tative results for a broad variety of cases suggests that 
some properties of many-body systems are founded not 
in detailed microphysics but on simple properties, specif- 
ically projection integrals and the relative dimensions of 
subspaces with good quantum numbers. 
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